W:  A  \o4  033. 


TECHNICAL 
.  LIBRARY 


TECHNICAL  REPORT  ARB RL-TR-02344 

CALCULATION  OF  BLOWOUT  GUN 
NOZZLE  TEMPERATURES 


Charles  S.  Smith 
J.  Richard  Ward 

July  1981 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 
BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


Approved  for  public  release;  distribution  unlimited. 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Secondary  distribution  of  this  report  by  originating 
or  sponsoring  activity  is  prohibited. 

Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.S.  Department  of  Commerce,  Springfield,  Virginia 

22161. 


The  findings  in  this  report  are  not  to  be  construed 
an  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 


The  use  of  tmde  names  or*  manufacturers  f  names  in  this  report 
does  not  constitute  indorsement  of  any  commercial  product . • 


_ 11NCLASSIE1ED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

TECHNICAL  REPORT  ARBRL-TR-02344 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

CALCULATION  OF  BLOWOUT  GUN  NOZZLE  TEMPERATURES 

5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

BRL  Technical  Report 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR^*; 

Charles  S.  Smith* 

J.  Richard  Ward 

8.  CONTRACT  OR  GRANT  NUMBERfaJ 

9-  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

H§  Army  Armament  Research  §  Development  Command 

US  Army  Ballistic  Research  Laboratory 

ATTN:  DRDAR-BLI 

.Aberdeen  Proving  Ground.  MD  2inns 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 

AREA  A  WORK  UNIT  NUMBERS 

1L162618AH80 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Armament  Research  §  Development  Command 

US  Ballistic  Research'  Laboratory 

ATTN:  DRDAR-BL 

Aberdeen  Proving  Ground.  MD  2 inns 

12.  REPORT  DATE 

JULY  1981 

13.  NUMBER  OF  PAGES 

46 

14.  MONITORING  AGENCY  NAME  A  ADDRESS^//  different  from  Controlling  Office) 

15.  SECURITY  CLASS,  (ol  thie  report) 

UNCLASSIFIED 

15a.  DECLASSI  FI  CATION/ DOWN  GRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thie  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  ebu  tract  entered  in  Block  20,  it  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


*Naval  Surface  Weapons  Center,  Dahlgren  Laboratory,  Dahlgren,  VA 

19.  KEY  WORDS  (Continue  on  reveree  aide  if  neceeemry  and  identify  by  block  number) 

Gun  Propellants  Calspan 

Gun  wear  and  erosion 
37mm  blowout  gun 
Bore  surface  temperature 
Nordheimfs  method 

ABSTRACT  (XZ antinua  mm  rmrerem  and  identify  by  block  number)  (  j  lc  ) 

Small  laboratory  devices  such  as  the  37mm  blowout  gun  at  the  BRL  have  a 
long  history  in  unraveling  the  factors  influencing  gun  barrel  wear.  Such  devices 
are  particularly  useful  when  testing  scarce  and  expensive  experimental  pro¬ 
pellants.  An  unresolved  question  is  whether  relative  propellant  erosivity 
measured  in  such  laboratory  devices  correlates  with  large-caliber  gun  wear. 

(Continued  on  next  pge) 

w>  A— »  W73  EO,„o»or,»ov..i,o«™..T«  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


_ UNC  LASS  TFT  FT) _ 

SECURITY  CLASSIFICATION  OF  THU  RAOE(TWli  Datm  Knitted) 


20.  ABSTRACT  (ContdL 

To  try  to  answer  the  question,  Nordheim's  scheme  for  computing  heat 
transfer  in  guns  and  erosion  yents  was  used  to  calculate  bore  surface 
temperatures  in  a  17.3  mm  diameter  nozzle  for  a  series  of  five  propellants 
for  which  the  wear  had  already  been  measured.  Empirical  expressions  for 
estimating  wear  in  guns  suggested  the  wear  should  increase  exponentially 
with  peak  bore  surface  temperature  if  the  nozzles  mimicked  wear  in  guns. 

A  linear-least  squares  of  natural  log  of  wear  vs .  peak  temperature  showed 
the  wear  in  nozzles  could  be  fit  to  such  an  expression,  specifically 
dn  a)  =  -  a.7  +  0.0094  T,  where  T  is  the  peak  surface  temperature  in  Kelvin 
and  wear  is  expressed  as  y/shot. 

To  test  how  well  the  dependence  of  wear  on  bore  surface  temperature  com¬ 
pares  to  large  caliber  guns,  bore  surface  temperatures  computed  by  Calspan 
from  heat  input  measurements  in  the  155m  Ml 85  cannon  were  plotted  against 
wear.  The  slope  was  0.0076  K“1  suggesting  the  blowout  gun  wear  may  exagger¬ 
ate  relative  wear  expected  in  a  large-caliber  gun  by  ten  percent. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEftFh.n  Dmf  Enff*) 


TABLE  OF  CONTENTS 


Page 

LIST  OF  TABLES . 5 

I.  INTRODUCTION  .  7 

II.  INPUT  PARAMETERS  AND  COMPUTER  CODES  .  8 

III.  RESULTS  AND  DISCUSSION . 10 

IV.  CONCLUSIONS . 17 

REFERENCES  . 21 

APPENDIX . 23 

DISTRIBUTION  LIST . 43 


3 


LIST  OF  TABLES 


Page 

1.  Compositions  and  Grain  Dimensions  of  M30  and  HFP  Propellants.  9 

2.  Compositions  and  Grain  Dimensions  of  M5,  M8 ,  and  Ml  Propel¬ 
lants  . 9 

3.  Propellant  Parameters . . 11 

4.  Experimental  Data  From  Blowout  Gun  Experiments . 12 

5.  Computed  Heat  Transfer  Results . 13 

6.  Comparison  Between  Experimental  Wear  and  Wear  Computed  With 

Equation  (5) . 14 

7.  Bore  Surface  Temperatures  and  Wear  for  155mm  Charges  Without 

Additive . . 15 

8.  Calculated  Heat  Input  and  Surface  Temperatures  for  3"/70  and 

37mm  Blowout  Gun . . . . . 17 


5 


I.  INTRODUCTION 


Small  laboratory  devices  such  as  vented  chambers^or  blowout  guns 
have  a  long  history  in  understanding  gun  barrel  wear.  Compared  with 
large-caliber  gun  tests,  such  devices  need  less  propellant,  need  less 
room  for  safe  firing,  and  have  inserts  which  can  be  removed  for  post¬ 
firing  analysis.  An  unanswered  question  is  how  to  relate  wear  in  the 
devices  with  wear  in  large-caliber  guns.  A  recent  .JANNAF  workshop  con¬ 
cluded  that  laboratory  devices  cannot  evaluate  wear-reducing  additives, 2 
and  that  test  devices  in  three  laboratories  have  given  contradictory 
results  with  the  same  propellants. 3 >4 

This  report  attempts  to  relate  erosivity  from  propellants  in  the 
BRL  37mm  blowout  gun  to  wear  expected  in  large  caliber  guns.  Empirical 
formulas^  6  presume  the  wear  is  exponentially  proportional  to  the  maxi¬ 
mum  bore  surface  temperature.  Bore  surface  temperatures  in  the  empira- 
cal  models  are  estimated  with  interior  ballistic  parameters  such  as 
charge  weight,  flame  temperature,  gun  diameter,  and  peak  pressure. ' 


1 

"Hypervelocity  Guns  and  the  Control  of  Gun  Erosion, "  Summary  Technical  ■ 
Report  of  Division  1 ,  NDRC,  Volume  1,  Washington,  DC  1946. 

2 

J.A.  Lannon  and  J.R.  Ward ,  "Workshop  Report  on  Wear-Reducing  Additives 
and  Their  Performance  in  Guns, "Proceedings  of  the  17th  JANNAF  Combustion 
Meeting,  CPIA  Publication  329 ,  November  1980. 

3 

A.J.  Bracuti,  L.  Bottei,  J.A.  Lannon,  and  L.H.  Caveny,  "Evaluation  of 
Propellant  Erosivity  with  Vented  Chamber  Apparatus,  "  1980  JANNAF  Pro¬ 
pulsion  Meeting,  CPIA  Publication,  315,  March  1980. 

^ J.R.  Ward,  R.W.  Geene,  A.  Niiler,  A.  Rye,  and  B.B.  Gr oilman,  Blowout 
Gun  Erosivity  Experiments  with  Double-Base,  Triple-Base,  and  Nitramine 
Propellants,  1980  JANNAF  Combustion  Meeting  CPIA  Publication ,  315 , 

March  1980. 

5j. M.  Frankie  and  L.R.  Kruse,  "A  Method  for  Estimating  the  Service  Life 

of  a  Gun  or  Howitzer,"  BRL  Memorandum  Report  No.  1852,  June  1967. 

6  (AD  #818348)  . 

C.S.  Smith  and  J.S.  O'Brasky,  "A  Procedure  for  Gun  Barrel  Life  Estima¬ 
tion,  "  Proceedings  of  the  Tri-Service  Gun  Tube  Wear  and  Erosion  Sympo¬ 
sium,  ADPA ,  Dover ,  NJ,  March  1977. 

7 

J.  Comer,  Theory  of  the  Interior  Ballistics^  of  Gunsj  John  Wiley  &  Sons , 
Inc.,  NY,  1950. 

Q 

L.W.  Nordheim,  H.  Soodak,  and  G.  Nordheim,  "Thermal  Effects  of  Propellant 
Gases  in  Erosion  Vents  and  Guns,"  NDRC  Armor  and  Ordnance  Report  No. 
A-262,  March  1944. 
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This  report  uses  Nordheim’s  scheme  to  compute  the  peak  surface 
temperatures  during  wear  tests  of  five  propellants.9 10  If  the  relative 
wear  of  the  propellants  in  the  blowout  gun  is  the  same  as  in  large  guns, 
then  a  plot  of  natural  logarithm  of  wear  v£  peak  temperature  should  be 
linear,  with  the  slope  identical  to  that  in  guns.  Wear  tests  ip  Reference 
9  incorporated  Niiler's  radioactive  technique^-9  to  correlate  mass  loss  and 
diameter  change. 

II.  INPUT  PARAMETERS  AND  COMPUTER  CODES 

Three  computer  codes  were  used  to  perform  the  calculations  dis¬ 
cussed  in  the  report.  The  first  program,  IB  3/70,  computes  interior 
ballistics  data  for  guns.  The  second  program,  NOZZLE,  computes  interior 
ballistic  data  for  a  blowout  gun  with  an  erosion  nozzle.  The  third  pro¬ 
gram,  NEWNSN,  computes  heat  transfer  using  files  of  IB  3/70  and  NOZZLE. 
Listings  for  each  program  are  contained  in  the  Appendix. 

The  assumed  heat  input,  as  described  fully  in  Reference  8,  was 
taken  to  be 


Q  =  JgXCppUAT  ,  (1) 

2 

where  q  =  heat  flux,  J/m  *s, 

X  =  friction  factor, 

C  =  specific  heat,  J/kg*K, 

3 

p  =  propellant  gas  density,  kg/m  3 
U  =  gas  velocity,  m/s,  and 

AT  =  temperature  difference  between  propellant  gas  and  wall,  K. 

Following  Nordheim's  method,  a  friction  factor  of  1/253  was  used 
for  the  17.3  mm  diameter  nozzle. 

Computations  were  done  for  the  five  propellants  evaluated  in  Re¬ 
ference  9.  The  compositions  are  listed  in  Tables  1  and  2. 


9  R.W.  Geene ,  J.R.  Ward 3  T.L.  Brosseau ,  A.  Wilier ,  R.  Birhnire,  J.J. 
Rocchio,  "Erosivity  of  a  Nitramine  Propellant , "  BRL  Technical  Report 
No.  020943  August  1978.  (AD  #A060590) 

10S.E.  Caldwell  and  A.  Niileri  "The  Measurement  of  Wear  from  Steel  Using 
the  Radioactive  bbCo3"  BRL  Report  No.  1923  3  September  1976.  (AD  #A030262) 
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TABLE  1.  COMPOSITIONS  AND  GRAIN  DIMENSIONS  OF  M30  AND  HFP  PROPELLANTS 


M30  HFP 

PPL-A-6372  PPL-A-6380 


Nitrocellulose  (12.6%N) 

28.0% 

29.3% 

Nitroglycerin 

22.5 

22.7 

Nitroguanidine 

47.7 

5.0 

RDX 

36.5 

Dioctylphthalate 

5.0 

Ethyl  Centralite 

1.5 

1.5 

Cryolite 

0.3 

Total  Volatiles  (Residual) 

0.2 

0.3 

Grain  Length,  mm 

7.78 

10.58 

Grain  Diameter,  mm 

1.59 

2.37 

Grain  Perf.  Diameter,  mm 

0.46 

0.77 

Grain  Web,  mm 

0.56 

0.80 

Grain  Geometry 

SP 

SP 

TABLE  2.  COMPOSITIONS  AND 

GRAIN  DIMENSIONS 

OF  M5 ,  M8 ,  AND  Ml 

PROPELLANTS 

M5_ 

M8 

Ml 

Nitrocellulose  (13.25%N) 

81.95% 

52.15% 

85.00% 

Nitroglycerin 

15.00 

43.00 

Ethyl  Centralite 

0.60 

0.60 

Barium  Nitrate 

1.40 

Potassium  Nitrate 

0.75 

1.25 

Diethylphtalate 

3.00 

Dinitrotoluene 

10.00 

Dibutylphthalate 

5.00 

Diphenylamine,  Added 

1.00 

Ethyl  Alcohol,  Residual 

2.30 

0.40 

0.75 

Water,  Residual 

0.70 

0.50 

Graphite 

0.30 

Grain  Length,  mm 

10.58 

25.4 

8.26 

Grain  Diameter,  mm 

3.92 

12.7 

3.68 

Grain  Perf.  Diameter,  mm 

0.41 

0.37 

Grain  Web,  mm 

0.69 

0.56 

0.64 

Grain  Geometry 

7  Perf 

Strip 

7  Perf 
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Propellant  parameters  required  for  the  computations  are  listed  in 
Table  3. 

The  term,  B/W,  is  used  to  relate  burning  rate,  r,  with  chamber  pressure,  P, 
and  propellant  web,  W,  by 

r  =  (P  +  44.8)  (|)  •  W  .  (2) 

Propellant  webs  are  listed  in  Tables  1  and  2. 

3 

For  the  blowout  gun,  chamber  volume  and  nozzle  area  were  315  cm 
and  2.36  cm2,  respectively. 


III.  RESULTS  AND  DISCUSSION 

The  experimental  results  from  Reference  9  are  listed  in  Table  4. 

Table  5  gives  calculated  time  to  maximum  pressure,  peak  nozzle  sur¬ 
face  temperature,  heat  input,  and  time  from  nozzle  rupture  to  peak  nozzle 
surface  temperature. 

A  linear  least  squares  fit  tested  the  exponential  dependence  of  wear 
with  peak  surface  temperature. 

£n  (w)  =  a  +  bT  ,  (3) 


where  w  =  wear,  y/shotj 

T  =  peak  bore  surface  temperature,  K,  and 
a,b  =  constants. 

The  best  fit  values  of  a  and  b  were  -9.7  ±  2  and  0.0094  ±  0.002, 
respectively,  with  error  representing  one  standard  deviation.  Table  6 
compares  experimental  wear  and  wear  computed  with  Eq.  (3).  Eq.  (3) 
clearly  underestimates  wear  for  Ml  with  a  peak  temperature  of  773  K. 
Frankle-Kruse  also  noted  their  expression  was  invalid  below  900  K; 
Smith-0' Brasky  report  a  similar  threshold  of  750  K. 

A  check  on  how  wear  varies  with  peak  bore  surface  temperatures  in 
guns  can  be  made  with  bore  surface  temperature  computations  by  Vassallo 
and  coworkers  at  Cal span  Corporation! 1  in  the  M185  cannon.  The  Cal span 
workers  compute  peak  bore  surface  temperatures  from  measured  total  heat 
input. 12  Recent  work  in  a  shock-tube  gun  with  different  metals  showed 
temperature  calculations  starting  with  measured  total  heat  input  could 
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TABLE  3.  PROPELLANT  PARAMETERS* 
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TABLE  4.  EXPERIMENTAL  DATA  FROM  BLOWOUT  GUN  EXPERIMENTS 
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TABLE  5.  COMPUTED  HEAT  TRANSFER  RESULTS 
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TABLE  6 


COMPARISON  BETWEEN  EXPERIMENTAL  WEAR  AND  WEAR  COMPUTED 
WITH  EQUATION  (3) 


Propellant  Peak  Temp. ,  K  Wear  exp't,  y/Shot  Wear  Cal’d,  n/Shot 


Ml 

Ml 


773 

875 


0.2 

0.1 


0.0005 

0.2 


M30 

966 

0.5 

M30 

1,130 

0.6 

M30 

1,228 

3.8 

0.5 

2.5 

6.2 


HFP 

968 

0.5 

HFP 

1,150 

1.1 

HFP 

1,235 

6.8 

0.5 

3.0 

6.6 


M5 

982 

0,8 

M5 

1,155 

4.1 

M5 

1,250 

18.5 

0.6 

3.1 

8.3 


M8 

1,107 

M8 

1,210 

MS 

1,407 

2.8 

9.6 

48.6 


2.0 

5.0 

32.1 
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predict  correctly  whether  the  melting  point  of  the  metal  was  exceeded. 
Starting  conditions  in  the  shock-tube  gun  were  varied  to  produce  tempera¬ 
ture  bracketing  the  metals'  melting  temperatures. 

Table  7  shows  bore  surface  temperatures  computed  for  three  155mm 
propelling  charges  which  have  no  wear-reducing  additive.  Wear  data  is 
available  for  the  M119  and  the  XM201E2  charges, 14 >15  Wear  sensor  data  in 
the  M185  cannon  showed  the  XM208  minus  its  additiye  was  five  times  as  ero¬ 
sive  as  the  XM201E2.  Wear  of  the  XM208  was  presumed  to  be  five  times 
greater  than  the  measured  wear  of  the  XM201E2  charge.  A  least-squares 
fit  of  the  three  points  in  Table  8  to  Eq.  (3)  produced  values  of  a  and  b 
of  -7.4  and  0.0076K-!,  respectively. 

TABLE  7.  BORE  SURFACE  TEMPERATURES  AND  WEAR  FOR  155mm  CHARGES  WITHOUT 
ADDITIVE 


Charge  Peak  Temp,  K  Wear,  p/shot  Ref. 


M119 

961 

0,9 

15 

XM201E2 

1,100 

2.6 

16 

XM208  (no  additive) 

1,306 

12.5 

12 

Figure  1  plots  wear  vs  peak  bore  temperature  for  the  nozzle  and  the 
155mm  gun.  Figure  1  suggests  relative  wear  among  propellants  will  be 
overestimated  compared  with  wear  in  guns  by  yen  percent.  Nonetheless, 
the  agreement  seems  close  enough  to  suggest  relative  propellant  erosivity 
can  be  inferred  in  nozzles. 


11 

F.A.  Vassallo,  "An  Evaluation  of  Heat  Transfer  and  Erosion  in  the  155 
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WEAR,  fx  /SHOT 
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5.0 


1.0 
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0.8  0.9  1.0  1.1  1.2  1.3 

PEAK  BORE  SURFACE  TEMP,  K,  *103 

Figure  1.  Wear  vs_  Peak  Bore  Surface  Temperature 
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Some  other  observations  from  the  calculations  of  bore  surface  tem¬ 
peratures  can  be  made  from  comparisons  of  bore  surface  temperature  vs. 
time  for  the  37mm  blowout  gun  and  a  3"/70  Naval  gun  depicted  in  Figures 
2-4.  Table  8  summarizes  pertinent  heat  transfer  data  for  the  three 
figures . 


TABLE  8.  CALCULATED  HEAT  INPUT  AND  SURFACE  TEMPERATURES  FOR  3"/70 
AND  37mm  BLOWOUT  GUNS 


Gun 

Nozzle 
Diameter , mm 

1/A 

Propellant 

Charge 

Mass,g 

Flame 

Temp,K 

Surface 
Temp , K 

Heat  In-^ 
put, J/mm 

37mm 

17.3 

253 

HFP 

89.8 

3,255 

1,235 

0.78 

3"/70 

- 

280 

Picrite 

3991. 

2,065 

1,224 

1.13 

37mm 

13.3 

239 

Ml 

159 

2,417 

1,254 

1.21 

The  first  point  of  interest  is  the  short  heating  time  of  the  blow¬ 
out  gun  in  Figure  2  vs.  the  3"/70  Naval  gun.  As  a  result  of  the  longer 
heating  time,  the  3"/70  gun  reaches  the  same  peak  bore  surface  tempera¬ 
ture  with  a  much  cooler  propellant.  The  effect  of  the  longer  heating 
time  is  also  reflected  in  the  larger  total  heat  input  for  the  gun.  The 
figures  also  point  out  why  propellants  with  flame  temperatures  near 
3800K  would  be  needed  to  get  the  wear/shot  characteristic  of  tank  guns 
with  the  17.3  mm  nozzle.  Figure  4  shows  that  Ml  propellant  produces  an 
equivalent  peak  temperature  as  HFP  if  charge  mass  is  increased  and  nozzle 
diameter  is  decreased.  A  new  chamber  would  be  required  in  the  37mm  gun 
to  be  able  to  load  160g  of  propellant. 


IV.  CONCLUSIONS 

1.  Bore  surface  temperature  vs.  time  were  computed  for  17.3  mm  diameter 
nozzles  in  the  BRL  37mm  blowout  gun  in  order  to  test  dependence  of  wear 
vs.  peak  bore  surface  temperature. 

2.  The  natural  logarithm  of  wear  from  the  nozzle  was  shown  to  be  linearly 
dependent  on  peak  bore  surface  temperature  which  agrees  with  empirical 
wear  models  and  with  calculations  of  peak  bore  surface  temperatures  by 
Calspan  Corp.  for  the  155mm  M185  cannon.  It  was  shown  the  slope  is  10% 
steeper  in  the  nozzle  suggesting  relative  propellant  erosivity  with  the 
17.3mm  diameter  nozzle  may  be  exaggerated  slightly  compared  to  large 
caliber  guns. 
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Figure  2.  Temperature  vs.  Time  In  37mm  Blowout  Gun  (17.3mm 
Diameter  NozzTe-HFP  Propellant). 
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TEMPERATURE  (DEG.  K ) 

400  G00  900  1000  12O0  14D0 


Figure  3.  Temperature  vs.  Time  In  3"/70  Gun  With  Plcrite  Propellant. 
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TEMPERATURE  (DEG.  K) 

400  600  600  100Q  1200  1400 


Figure  4.  Temperature  vs.  Time  in  37mm  Blowout  Gun  (13.3mm 
Diameter  NozzTe-Ml  Propellant). 
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APPENPIX 


COMPUTER  CODRS 


COMPUTER  CODES 


The  first  program,  "IB3/70",  computes  interior  ballistics  data 
for  guns . 

The  second  program,  "NOZZLE",  computes  interior  ballistics  data 
for  a  blowout  gun  with  an  erosion  nozzle. 

The  third  program,  "newnsn"  computes  heat  transfer  and  conduction 
data.  The  program  "newnsn"  uses  information  stored  on  files  !'DATA" 
and  "DATA1"  by  "IB3/70"  or  "NOZZLE". 
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"IB3/70" 


Input  data  begins  in  line  2860.  Data  to  be  stored  in  each  line 


is  as  follows: 

Line 

Data 

2860 

Projectile  weight  (lbs) 

2870 

Projectile  travel  (inches),  chamber  volume  (cu.in.) 

2880 

Pidduek-Kent  constant,  integration  step  size  (sec) 

2890 

Estimated  projectile  velocity  (ft/sec),  bore  diameter (in) , 
bore  area  (sq.in.) 

2900 

Shot  start  pressure  (psi) 

2910 

Igniter  weight  (pounds) 

2920 

Igniter  impetus  (in  lbf/lbm),  ratio  of  specific  heats, 
flame  temperature  (°K),  covolume  (cu.in. /lb) 

2930 

Propellant  weight  (lbs) 

2940 

Propellant  force  (  in  lbf/lbm) ,  ratio  of  specific  heats, 
flame  temperature  (°K),  covolume  (cu.in. /lb),  density (lb/cu. in 

2950 

Propellant  grain  diameter  (in),  length  (in),  perforation  dia¬ 
meter  (in) ,  and  number  of  perfs 

2960 

Propellant  identification  (an  alphanumeric  string) 

2970 

PI,  B (where  the  burning  rate  (in/sec)  is  B* (space  mean 
pressure  in  psi  +  PI)),  number  of  integration  steps  per 
printout,  expected  peak  chamber  pressure  (psi),  expected 
peak  acceleration  (g's). 

2980 

Number  of  "resistive  pressures'1  listed,  followed  by  the 
appropriate  number  of  pairs  of  travel  distance (in)  and 
resistive  pressure  (lb/sq.in.) 

2990-3020 

These  lines  are  used  by  the  program  to  aid  in  calculating 
the  fraction  of  the  7  perf  grain  burnt  after  splintering. 

3025 

Use  this  line  to  provide  data  if  called  for  by  the  input 
statements 

Data  requested  by  input  statements: 

"Number  of  parameters  to  be  varied":  Enter  a  number  (0  to  5). 

"Number  of  proj  wt,  cham  vol,  pmax,  travel  distance,  charge  wt":  The 
number  of  non-zero  enteries  should  be  as  indicated  in  the  immediately 
previous  input.  The  values  for  the  parameters  indicated  should  have  been 
listed  in  line  3025  of  the  program. 

"Option  No  2" 

Enter  a  number  (0  to  7).  The  program  will  then  perform  as  indicated 
below: 
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NUMBER  ENTERED 


ACTION 


0  Compute  IB  trajector  for  input  data  (i.e.,  data 

in  lines  2860-2980) 

1  Change  B  to  match  indicated  peak  pressure 

2  Change  charge  weight  to  match  indicated  pressure 

3  Change  B  to  match  indicated  acceleration 

4  Change  charge  weight  to  match  indicated  acceleration 

5  Change  charge  weight  to  match  indicated  velocity 

6  Print  out  a  blurb  of  the  above 

7  Same  as  0,  but  program  writes  data  on  files  DATA 
DATA1  for  future  use. 

Note  that  if  the  first  input  requested  was  not  zero  multiple  runs  will  be 
performed.  This  is  not  useful  if  "Option  number  2"  is  7. 
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10. 

OP 'II  OH  ERSE  1 

20  O’ 

'ERL  HP 

30 

DIM  C < 50 > 

40 

DIM  P  *  [ 10] 

50 

DIM  T  <  2 , 1 0  > , F  <  2 , 1 

2  >  ,  R  *s  6 , 3  > 

60 

Eheat-O 

70 

LET  26=0 

SO 

D E F  F N fl <  Y 9 >  =  IN T  < Y 9*1 U 0 G > / 1 0 0 0 

90 

F:  E  A  D  G 1 ,  G  2 ,  G  3 ,  G  5 , 

G,  6 ,  G  7 ,  i  j  3 ,  G  4  , 

100 

LET  E7=G7*12 

110 

READ  I  1  , 12,  13, 14, 

15 

120 

READ  P 1 , P 2 , P 3 , P 4 , 

P5,P6,P7,P8, 

130 

READ  Ft 

140 

READ  H , B , M 2 , M 3 , A  9 

150 

IF  PO< >7  THEN  190 

160 

LET  W 1  =  < P 7 - 3 * P 9 > 

4 

17G 

LET  W2=  <  < P7-W1 >  A2 

■"  7  *  <  p  9  +  N 1  )  2 

ISO 

GOTO  200 

190 

LET  N 1  —  1 0 0 0 

2  0  0 

READ  K 

216  FOR  1=1  TO  K 

220  READ  TCi,I  >,  T.<2,  I  > 

236  NEXT  I 
246  FOR  1=1  TO  12 
25U  READ  F(1JI),F(2,D 
260  NEXT  I 

270  PRINT  " n UMBER  OF  PARAMETERS  TO  BE  VARIED" 

280  INPUT  D 

290  IF  D=0  THEN  516  * 

3O0  PRINT  "NO  OF  FROJ  NT  ,  CHAM  VOL  ,  PMRX ,  TRAVEL  III  ST,  CHARGE  NT" 

310  I NPUT  D 1 , D2 , D3 , D4 , D5 

320  IF  D1=G  THEN  366 

330  FOR  N 1 = 1  TO  D1 

346  READ  LOU). 

350  NEXT  N1 
360  IF  D 2  =  O  THEN  408 
378  F 0 R  N 2=1  TO  D  2 
336  READ  MCN2) 

390  NEXT  N2  ; 

400  IF  D3=0  THEN  470 
410  FOR  N3= 1  TO  D3 
420  READ  GKN3> 

430  NEXT  N3 

440  FOR  N 4  = 1  TO  D4 

450  READ  N < N4 ) 

460  NEXT  N4 
470  IF  D 5  =  0  THEN  510 
480  FOR  N 5  =  Q  TO  D5 
490  READ  Y < N 5 > 

5O0  NEXT  N5 

510  PRINT  "OPTION  NO  2  <OOX>-5>  ;  INPUT  6  fOR  EXLANAT  ION" 

520  INPUT  Cl 

530  IF  C 1 < 6  THEN  760 

540  IF  Cl =6  THEN  620 

550  ei=o 

560  CREATE  "DATA" ,126 
570  CREATE  "DATA1" ,3 

580  ASSIGN  #1  TO  "DATA" 

581  BUFFER  #1 

590  ASSIGN  #2  TO  "DATA1" 

60O  Eh eat = 1 
610  GOTO  930 
620  PRINT 
630  PRINT 
640  PRINT 

650  PRINT  "TO  MATCH"  ;  TAB <20 )  ;  "BY  ADJUSTING'1 ;  TRB<41J> :  11  INPUT" 
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660 

PRINT 

670 

PRINT  “NOTHING" ; TAB (4 

680 

PRINT  " PM AX" ; T ABC 20) ; 

690 

PRINT  " P M A X “ ; TAB C 20) ; 

700 

PRINT  “ACCEL" ; TAB <29 ) 

716 

PRINT  "ACCEL" ; TAB C 20) 

720 

PRINT  "VELOCITY" ; TAB  C 

730 

PRINT 

740 

PRINT 

750 

INPUT  Cl 

760 

IF  B=0  THEN  930 

770 

IF  Dl=8  THEN  800 

( y  0 

FOR  N 1 = 1  TO  D1 

790 

LET,  G 1  =  L  C  N 1  ) 

800 

IF  B 2  =  0  THEN  830 

8  1  O 

FOR  H2= 1  TO  B2 

820 

LET  G 3  =  M < N 2 > 

830 

IF  B3=0  THEN  860 

840 

FOR  N3= 1  TO  D3 

850 

LET  M 3  =  Q C H 3 > 

860 

IF  D 4  =  0  THEN  890 

870 

FOR  N4= 1  TO  B4 

8  y  o 

LET  G2=NCN4) 

890 

IF  D 5 = 0  THEN  930 

900 

FOR  N5  = 1  TO  B5 

910 

LET  P 1 = Y C N 5 ) 

9  2  0 

NEXT  N5 

9  3  0 

LET  J1 =1 1*12/(13-1 > 

940 

LET  J  2  =  J  1 / 1 4 

950 

LET  J 3  = I 1*12 f I 4 

9  6  O 

LET  J5=I1*I5 

970 

LET  Q 1 =  P 1  * F 2 s C P 3-1 ) 

986 

LET  Q 2  =  Q 1 / P 4 

9  9  0 

LET  Q3=P1*P2/P4 

1  000 

LET  Q 4  =  P 1  / F* 6 

1010 

LET  Q 5  =  P 1  * P 5 

1020 

LET  Q 6  = I 1 +  P 1 

1030 

Q7  =  P  8  *<P?  2  -  P  G  *  P  9  -  2  ) 

1040 

LET  Q 8  =  G 6 / 2 

s ) ;  “  0 " 

'BETA” ; TAB < 43); "l" 

‘CHARGE  WT";TAB<43>; "2“ 

;  ,,BETR,,;TfiB<43>;  "S" 

!  “CHARGE  WT,,;TAB<43>;  "4" 

>0) ; “ CHARGE NWT" ; TAB < 43) ; “5" 


105O  LET  CU  =  .  38*G8M .  5* < G2+G3/G4 ) *<P4-380 ) 

106G  LET  01=Ol/<  1  +  .  6*G8A2 .  175/PI*.  3375) /{.  1  44*G7*G7) 
1G70  LET  H 1 =0 1  +  CL  1 +  Q6/G5  > / 7 7 2 
1030  LET  0 2=1 / C P 3 - 1 ) 

1090  LET  H2=i+Q6/CG1*G5> 


1100  LET  O 3=1/ C < 2*02+3 )/G5+2*< 02+1 >*G1'Q6> 

1110  LET  H3=< 1-03>M-Q2~1 > 

1120  LET  H4  =  G4*386  yG 1 
1130  IF  M2 >=150  THEN  1160 

114G  PRINT  “TIME  BIST  PRES  VEL  TEMP 

1150  PRINT  "FR.  ACCEL  BASE" 

1160  LET  T=1000 
1170  LET  T2=T<2 , 1 ) 

1130  LET  H - 0 
1190  LET  M 1 = 1 

1200  LET  R C 1 , 1 ) =  R C 1 , 2 )  =  R < 1 , 3>=B0=Z-T2«S=B3=B4=B5=Z5*Z 
1210  LET  X1=RC1,2> 

1220  LET  V 1 =R C  1 , 3) 

1230  FOR  1=1  TO  3 
1240  LET  R  C  2 ,  I  )  =  R  C  1  ,  I > 

1250  NEXT  I 
1260  LET  J  =  2 
1278  GOSUB  2230 
1230  LET  S-S+Q8 
1290  LET  J  =  3 
1300  FOR  1=1  TO  3 
1310  LET  R<2, I)=R<2,I>*G6 
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1320  LET  R  •;  3  ,  I  )  -  R  <1,1  >  +  R  C  2 ,  I  >  /  2 

1336  NEXT  I 

1  3  4  0  G  0  S  U  B  2  2  8  8 

1350  LET  J  =  4 

1360  FOR  1=1  TO  3 

1370  LET  R ( 3 ,  I >  =  R  <  3 , I )  * G 6 

1380  LET  R  <  4 ,  I  >  * R  <  1  ?  I  >  +  R < 3 ,  I  )  /2 

1390  NEXT  I 

1  4  0  0  G  0  S  U  E  2 2  3  O 

1410  LET  J=5 

1420  LET  S=S+Q8 

1430  FOR  1=1  TO  3 

1440  LET  R<4, I>=R<4, I>*G6 

1450  LET  R <  5 .  I  )  =  R  <  1  ,  I  >  +  R  <  4  ,  I  ) 

1460  NEXT  I 

14  7 0  GOSU B  2280 

1480  FOR  1=1  TO  3 

1490  LET  R<5, I)=R<5, I )*G6 

1508  LET  R(l, I ) =R < 1 , X>  +  < R<2, I>+R<3, I>*2+2*R<4, I>>/6 
1510  LET  R<6, I )=R< 1 , I > 

1520  NEXT  I 
1530  LET  J=6 

1548  IF  Eheat  THEN  PRINT  8 i ; V, X , P 1 *?/< G3+G4*H-P 1  * U ~Z ) /R6> , J 
1550  IF  M2*INT<Ml/M2X>Mi  THEN  1590 

1 560  PR  1  NT  S* 1 800 ; TRB < 8 > ; FNR < R <1 , 2 ) > ; TftB < 1 8  > ; I NT  < BQ > ; TAB < 27 ) 5 
1 57G  PRINT  I  NT  < R U „ 3  > ✓ 1 2  > ; TRB < 34 > ; I  NT < T  ) ; TRB  C  45  > ; FHR< Z > 5 
1580  PRINT  TRB < 54); I  NT <R 1/386. 09> ; TRB < 63) ; INT<B2> 

1590  IF  R < 1 , 2 ) >  =  G 2 - 2 5  THEN  1610 

160O  GOTO  1630 

1610  H  =  N  + 1 

1626  LET  C < N >  =  B2 

1630  IF  R < 1 , 2 ) > = G 2  THEN  1730 

1640  LET  Mi =M1 +1 

165G  LET  V=R< 1,3) 

1668  IF  B 0 <  = B 3  THEN  1210 

1679  LET  fl 5  =  fl 1 / 3 8 6 . 0 9 

1680  LET  B3=B0 

169Q  LET  B 8  =  B 2  ; 

1700  LET  B 4  =  8 
1718  LET  E  5  =  R  <  i  ,  2 ) 

1720  GOTO  1218 

1730  LET  X=<G2-X1  )/(R(l  , 2) -XI ) 

1740  LET  V  =  V 1+X*<R(1 , 3 ) - V 1 ) 

1750  LET  S 1 =  S - Q 8 

1760  LET  29  =  C < N ~ 1 ) +X*K C < N ) - C < N - 1  >  } 

1770  IF  C 1=0  THEN  I960 
1  7  8 0  B r  an c  h = C 1 

1798  IF  <Branch>-l )  RNB  <Branch<6)  THEN  ON  Branch  GOTO  1300, 1830, i960, 1890, 1920 
1800  IF  RES < B3-M3 ) <200  THEN  I960 
1810  LET  B= B /  < 1 - . 8  5*  <  M 3 - B 3  >  /  <.  M3 + B 3  >  ) 

1820  GOTO  1940 

1830  IF  RBS < B3-M3 ) < 200  THEN  I960 
1840  LET  P 1 =P i / < 1 ~ . 9*  < M3-B3 > / < M3+B3 ) > 

1850  GOTO  1948 

I860  IF  RBS<R9-R5><100  THEN  1960 
1870  LET  B  =  E / < 1-1 .  1  * < R9-R5 > s < R9  +  R5 ) ) 

1880  GOTO  1940 

1890  IF  RBS<R9-R5)< 100  THEN  I960 
19G0  LET  P  1  =P  1  /  1  -  .  9*  <  R9-R5  )  /  < R9  +  R5 )  ) 

1910  GOTO  1940 

1920  IF  RBS  (.  G7- V/Q2 )  <10  THEN  I960 

1938  LET  Pl=Pl/a-l  .5*(G7-V/12>/(G7+V/12>> 

1940  LET  E7=V 
1950  GOTO  930 
I960  PRINT 

1970  IF  Eheat,  =0  THEN  2000  20 


1980 
2000 
2010 
2O20 
203O 
2040 
2O50 
2060 
207O 
2080 
2090 
2100 
2110 
2120 
2130 
2 1  40 
2150 
2160 
2170 
2180 
2  1  90 
220O 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2288 
2290 

2309 

2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
248G 
24  90 
2500 
2510 
2520 
2530 
2540 
2550 
2568 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 


CALL  F-: j  <#1  ,  #2, P1/<G3+G2*G4>  ,  <V>  ,  S,G6,T  ,  P5,  G2+G3/G4 ,  G8 ,  G3.'G4  > 

IF  M2 >=50O  THEN  2190 

PRINT  " PRO J . NT .  = " ; G 1 ; "LBS . " 

PR  I  NT  “  I  HI  T  .  CHAM .  VOL .  =  " ;  G3 ;  "  CU  -  I N .  " 

PRINT  "CHG  NT  ; PI ; “LBS. " 

PR  I  NT  “TRAVEL  D I  ST .  =  "  ; G2 ; M I N . H 

PRINT  "MUZZLE  VELOCITY  = ,! ;  I  NT  <  V/ 1 2  > ;  "  F*'S 11 

PRINT  “MAX.  CHAM.  PrES.  = 11 ;  I  NT  <  E3  >  ;  "  PS  I  AT"  ;  £4;  "SEC.  OR"  ;  B5;  "  JN,  " 
PRINT  "MAX.  BASE  PRES .  =  " ; I NT <B8) ; "PS I " 

PRINT  "MUZZLE  PRESSURE  =  " ; I  NT < Z9 > ; "PS I  AT  H: SI; “SEC" 

PRINT  "MAX.  ACCEL  = " ; I  NT  < A5  > ; " GS " 

IF  Z<1  THEN  2136 

PRINT  "FRACTION  BURNT  =  "  ;  2 ;  "  AT 11 ;  Z5 ;  "IN  FROM  MUZZLE" 

GOT G-  2140 

PRINT  "FRACTION  BURNT  - " ; Z 

PRINT  "BETH" ; B 

PRINT 

PRINT  G9  s I  1 , T < 1 , 1 > , T < 2 , 1 > , T < 1 , 2) , T < 2, 2 > 

PRINT 

GOTO  2740 

IF  E9= 1  THEN  2240 

PRINT  " PRO J  NT  BIST  CHAM  VOL  CW  PMHXP 

PRINT  "IV  ACCEL  BETA" 

LET  E9= 1 
PRINT 

PR  I  NT  TAB ( 1 > ; G 1 ; TAB  < 8 > ; G2 ; TAB < 1 5 >  J G3 ; TAB < 24) ; P 1 ; T HB  <  3?  > : 

PR  I  NT  I  NT < B3 > ; T ; B <40 > ; I  NT < V/ 1 2 > ; TAB  <48 > ; I NT <  A5  > ; TAB <  56 > ; B 

PRINT 

GOTO  2740 

LET  U=R < J , 1 >  • 

LET  X  =  R < J , 2 > 

LET  V  =  R  <  J , 3 ) 

IF  2>=1  THEN  2350 
IF  U >  =  W 1  THEN  2400 

LET  Z  =  1  - <  P8-U  *  <  <  P7-U  >  A2-P0*  < F9  +  U )  -*  2  >  /Q7 
GOTO  2470 
LET  2=1 

IF  Z6= 1  THEN  2470  j 

LET  2  5 = G  2  - X 

LET  26=1 

GOTO  2470 

FOR  K = 1  TO  12 

IF  F  < 1 , K > * W 1 > U  THEN  2440 

NEXT  K 

GOTO  2350 

LET  U3=<  IJ/W 1  -  F  <  1  ,  K~1  >  > ✓  <  F  <1  ,  K >~F<  1  ,  K-  1  >  > 

LET  M4=  < P8-U > /P8 

LET  2= 1 - W 4 * <  F  < 2 , K - 1 >  +  <  F  <  2 , K > ~F  <  2  s  K-*  1 > )*W3>*W2 
LET  T 1 =  0 

IF  X  =  0  THEN  2560 

FOR  K=2  TO  10 

IF  X  <  =  T  <  1 ,  K  >  THEN  2530 

LET  T 1 =  T 1  +  < T  < 1 , K ) -T < 1 , K- 1  >  >* <T < 2 , K>  +  T < 2, K- 1 >  > /2 
NEXT  K 

LET  T2=T<2,K-1 )+a<2,K)-T<2,K-l))*(T(l,K)-X)/(T<l,K>-T<l 
LET  T  1  =  T  1  +  <  T  2  +  T  <  2 , K~ 1 > >  * < X-T < 1 , K- 1 > > /2 
LET  T 1 =T 1 *G4 

LET  T  =  <  J  1  +  Q 1  *  2  -  H 1  * V * V - T 1 >,'< J2  +  Q2*Z> 

LET  V0  =  G3  +  G4*X-Q4#  < i  -2  > - J5-Q5*Z 
IF  V 0 < 0  THEN  2608 
GOTO  2620 

PRINT  G3,G4,Q4, J5,Q5,Z 
GOTO  3030 

LET  B1=T*(Q3*2+J3)/V0 
LET  B 2  =  B 1  H 2 

LET  E0=B2*H3  30 
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2656 
266  8 
2670 
2680 
2690 
2708 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2868 
2370 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2  9  8  0 
2990 
3000 
3010 
3020 
3030 
3990 
4000 
4005 
4020 
406O 
4070 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
426u 
4  270 

4280 

4281 
4290 

4300 

4301 


IF  X<>8  THEN  2670 
IF  B1CG9  THEN  2690 
LEI  A1=H4*<B2-T2> 

GOTO  2700 
LET  A  1=8 

LET  R<J, 1>=2*B*<B1+A> 

LET  R < J , 2)= V 
LET  R<Jf3>=ftl 
RETURN 

IF  H=0  THEN  2850 
IF  D5=0  THEN  2770 
NEXT  N5 

IF  Ii 4 <  =  0  THEN  2790 
NEXT  N4 

IF  D3  =  0  THEN  2819 
NEXT  N3 

IF  D2=e  THEN  2830 
NEXT  N2 

IF  D 1  =  0  THEN  2850 

NEXT  HI 

PAUSE 

DAT A  15 

DATA  185.84,360 

DATA  3,. 00005 

DATA  3400, 3,7.21 

DATA  2500 

DATA  .02 

DATA  1 . 152E+06, 1 . 25,2000,30 
DATA  8.8 

DATA  3445900, 1 . 269,2065,32.81 , . 056- 
DATA  . 1425, . 35, .9132,7 
DATA  " EX 6586" 

DATA  6500, . 0O0856673, 280,52388,0 
DATA  2,8,1800,508,1800 

DATA  1,1,1 .086, 8. 967,  !  .  024, 8. 874,  1 . 05.5, 8.735 
DATA  1.101,0. 569, 1 . 162,0.395, 1 . 24,8.233 
DATA  1 . 292,0. 191,1. 347,8. 104,1. 489, 0,049 
DATA  1 . 48,8.814, 1 . 557,8 
STOP 

SUE  E  j  <  #  1 ,  #2 ,  Rho  i  ,  V  i  ,  T  i  ,  Dt. ,  T  e r»i p  i  ,  Ada,  Xm .  D  ,  C 1  > 

OPTION  EASE  1 

OVERLAP 

Eoa=Eob=0 

X  t  =Xm 

Kount =  T i /Dt 
V i = V i /Xm 
Rh= 1 /Rho i 
T  p  =  8 
G  m  1  = .  2  5 
C  u = . 34 

I  N P U T  "  C v ,  fid  a ,  G m  1  "  ,  C M ,  Ada,  Gro  1 
A 110: K o u  n t = K o u n t  + 1 
T  p  =  T p  +  Dt 
T=T i +Tp 
D s- n*i  =  1  +  Vi  *Tp 

V  =  V  i  /  Ii  e  m 
Rho = Rho i /Dem 

Temp=T  emp i  <Rh~Ada>/< 1 /Rho- Ada) ) A  Q m 1 

V  =  V*Xrn 

PRINT  #1 ; V, Xt , Rho, Temp 

IF  Kount  MOD  408  THEN  A118 

IF  Eoa  THEN  A40 

IF  V*Rho >  = 1 80  THEN  A110 

Eoa=  1 

Dt=4*D-t 

PR IHT  #1 ; 1E51 ,0,0,0  31 


CO  CO  CO 


4305  GOTO  ft  110 

4310  ft40:IF  Rho>=. 0O050  THEN  ft  1 10 
4315  PRINT  #1 ; 1E51 , 0, 0, O 
4320  X irt 
4330  Dt=Dt--4 

40  PRINT  #2;  Rhoi  ,Vi  ,Ti  ,-B* ,  T*tnp4  ,  flda,  Xm,  D,  C 1 
60  SUBEND 
70  END 


t 
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"NOZZLE" 


Input  data 
of  lines  20  and 

Parameter 

F0 

Beta,  PI 

C 

Rhop 

Ada 

TO 

Gamma 

Force 

C0,C2,C2P 

VO 

Ast 

Pstart 

G6 


is  lines  110  and  120  corresponds  to  the  READ  statements 
30.  Input  dat  is  as  follows: 

Data 

Initial  web  fraction  remainging.  (0<F0<1) 

If  F0  is  too  small  the  initial  pressure  will 
exceed  the  intended  nozzle  start  pressure. 

The  assumed  burning  rate,  in  inches  per  second, 
is  Beta  (P+P1)/W  where  P  is  the  chamber  pressure 
in  psi  and  W  is  the  web  size  (inches) 

Charge  weight  (pounds) 

Propellant  density  (lb/cu.in.) 

Propellant  covolume  (cu. in/lb.) 

Adiabatic  flame  temperature  (°K) 

Ratio  of  specific  heats 

Impetus  (in  lbf/lbm) 

These  are  the  kg,  k2,  and  -V-2  of  NDRC  A142  of 
5  Feb  1943;  they  are  used  to  help  calculate  the 
amount  of  propellant  that  has  been  burnt. 

Chamber  volume  (cu.in.) 

Area  of  the  erosion  nozzle  (sq.in.) 

Nozzle  start  pressure  (psi) 

Integration  step  size  (sec) 
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10  OVERLAP 

2  0  R  E  A  D  F  0 ,  B  e  ta,F'l,C,  R  h  op,  Ada,  T  0 ,  G  am  rn  a ,  F  o  r  ce,  CO,  C  2 ,  C  2  p 

SO  READ  V0 , Ast , Pst  art , G6 
40  S i =S=0 

50  Zeta=<2^<Gamrna+l  )  <  1  /(Garnma-  1 )  ) 

Go  Ep=<  1  -  2  *  Z  e  t  aO'Gamma 

?0  P h e- 1  p  =  <  2 2  <  G  am m a+ 1  >  >  A  <  G  am m a--  <  G  am  m a- 1 > ) 

80  0  8  -  G  6  -■  2 
90  Break  =0 

1 08  I  MAGE  BB . DDDD , DD . DDL , DDDD . DDUD , DDDDDD , BDDDDDDDD ,  DDL .  DDDDD ,  8D 
j 1 0  IlRT fl  ®  •  500 »  ■  005280 , 6500 ,  .  300 ,  .  057 ,  32 .81, 2065 ,  1 . 269 , 3445000 ,  .  8502 

120  DATA  19. 2, . 098, 41000, . 000050 

130  R  <  1  1  )  =  F  O 

140  R<1,2>=0 

150  Fount  =  0 

160  CREATE  “DATA", 126 

170  CREATE  " BA  T  A 1  ",  3 

180  ASSIGN  #1  TO  "DATA" 

190  ASSIGN  #2  TO  "BATR1 “ 

200  BUFFER  #1 

210  FOR  1=1  TO  2 

220  LET  R<2, I)=R< 1,1) 

230  NEXT  I 

240  LET  J=2 

250  GOSUB  660 

260  LET  S = S  +  Q 8 

270  LET  J=3 

280  FOR  1=1  TO  2 

290  LET  R<2, I >=R<2, I )*G6 

300  LET  R(3, I )=R(1 , I )+R<2, I 02 

310  NEXT  I 

320  GOSUB  660 

330  LET  J  =  4 

340  FOR  1=1  TO  2 

350  LET  R<3, I )=R<3, I )*G 6 

360  LET  R<4, I )=R<1 , I )+R<3, I 02 

370  NEXT  I 

380  GOSUB  660  •  1 

390  LET  J  =  5 

400  LET  S=S+Q8 

410  FOR  1=1  TO  2 

420  LET  R < 4 , I  ) =R  <  4 , I)*G6 

430  LET  R < 5 , I ) =  R < 1 , I  )  +  R  (  4 ,  I ) 

440  NEXT  I 

450  GOSUB  660 

460  FOR  1=1  TO  2 

470  LET  R(5, I )=R<5, I  )*G6 

480  LET  R  <1 ,  I  >  =R  <  1  ,  I  )  h-  <  R  <  2 ,  I  )  +R  <  3 ,  I  )  *2+2*R<  4  ,  I )  +R  < 5 ,  I  >  )  /6 
490  LET  R  <  6 ,  I >=RC 1 ,  I) 

500  NEXT  I 

510  PRINT  USING  100; S, 2, Si , T, P, Rho, Ust 

520  IF  U s X  =  0  THEN  210 

5  3  0  K  o  u  n  t,  =  K  o  u  n  t  + 1 

540  PR  I  NT  #  1  ;  List ,  1 ,  Rhost ,  Tst 

545  IF  Fount  MOD  400  THEN  210 

550  IF  NOT  Break  OR  < Rho >.8035)  THEN  210 

570  IF  SI  THEN  SI  a 

580  Sl=l 

590  PRINT  #1; 1E51, 0,8,0 
60O  G6=4"*G6 

605  08=4*08 

610  GOTO  210 

615  SI  a:  IF  Rho >.001  THEN  210 

620  PRINT  #1; 1E51, 0,0,0 

630  PRINT  #  2  J  0 ,  0 , 0 ,  G  6  4  ,  O ,  A  d  a ,  1,  1 . 25*SQR<  Ast  )  ,  1 
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1329, .367 
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640  L 0 fl D  " n e w n sn“ 

650  STOP 
660  F=R<J,1) 

670  Si=R<J?2> 

680  2=< l-F>*<C0+C2p*F> 


680  IF  2 > 1  THEM  2=1 

70O  IF  F >0  THEN  Z=< 1-F>*<C0tC2*F> 

710  IF  F<-.5  THEN  2=1 
720  Zhe  1  p=V0-C*  <  1  -2  >  ^Rhop 
780  Vpres*Zhelp-<C#Z~Si  !* * fl d a 
749  Vth=Zhelp-C*Z*flda 

7  bo  T=T0  *  (.  Vt>  h*  (  Z^S  i  /C  )  /(,  Vpres*Z)  )  6  ani  m  a  - 1  > 

760  P  =  F o r- c e * < C * 2 - S i  >*V<T0*Vpres> 

770  Break  = Break  OR  <P>Pstart> 

780  Rho=<C*Z-Si )/2help 


7  9  o  s  i  g  Hi  a=  R  h  o  /  (,  1  —  fl  d  a*  R  h  o  !j 


800  F; <  J ,  l  > =-Bet  a* < P+P 1 ) 

810  IF  NOT  Break  THEN  RETURN 

828  He  1 per  =  Ep*flda*Si gwa 

830  S  i  gmaf  t=2et-  a*S  i  gma*  C 1  +  He  1  per > 

840  RKost=Si  grnast./<  l+flda*Si  grnast ) 

8  5  o  T  s  t  =  2  e  T  *  <  1  +  (  6  am  m  a- 1  >  *  H  e  1  p  e  r  >  /  C  G  an)  w  a+ 1 ) 

8t>0  Pet  =Phe  1  p*P*  C 1  +GarHfiiaeHe  1  pep ) 

87o  Ust=8QR(386*6aft'irna*F,st*(  l+Rda*Rhost  >  ■■-S^Rhoet  ) 
880  R( J, 2>»flst*Uftt*Rhost 
890  RETURN 
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"newnsn" 


The  data  listed  in  this  program  should  not  normally  be  changed. 
Input  parameters  are  requested  as  indicated  below.  Default  values 
are  as  indicated. 


Message 


Keyboard  Entry  . 


0  FOR  GRAPHS,  1  FOR  NO  GRAPHS 
TIMM,  TEMPM 

ENTER  LABEL 
GUN  TEMPERATURE? 

ALCP 

XINT,  TINT 


NRDS,  NLA,  IPROPT,  LHELP 


INSIDE  DIAMETER, 
OUTSIDE  DIAMETER 


0  or  1 

Maximum  time  (sec)  and  temperature 
(°K)  to  be  shown  on  graph 
An  alphanumeric  label  for  the  graph 
Temperature  of  gun  steel 
(default  is  300° K) 

(default  is  0.4/(14.8  +  4*logio  D)^, 

D  is  the  gun  caliber  in  centimeters) 
Distance  from  idealized  breech  (inches) 
for  these  calculations,  time  to  stop 
calculations  or  go  to  next  round  (sec) 
(default  =  projectile  base  or  nozzle 
throat,  D/100) 

Number  of  rounds  in  this  burst,  number 
of  mathematical  layers,  print  option 
(0  for  full  print,  1  for  1/4  as  much, 

2  for  very  little  printing) ,  mathematical 
layer  to  use  for  bore  surface  between 
rounds  (default:  1,  7,  0,  4) 

Inside  and  outside  diameter  of  the 
gun  (inches) 
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1  OPTION  BASF  1 

5  DIM  Label $C 40] 

6  INPUT  "0  FOR  GRAPHS, ]  FOR  HO  GRAPHS ", Enq 

7  IF  Eng  THEN  355 

10  PLOTTER  IS  " GRAPHICS “ 

11  PEN  2 

12  DEG 

2  9  I  N  PUT  "TI N  M ,  T  E  M  P  M  " ,  T  i  m  m ,  T  *  m  p  m 

25  L INPUT  "ENTER  LABEL"  ,  Labe]  $ 

26  Label  =LEN<L*bel $> 

30  LOCATE 

4  0  SCALE  -  T  i  m  m  /  1 0 ,  T  i  m  rn  ,  2  0  O ,  T  ■=  m  p  m  +10  0 
45  CL  IR.  -T  i  mm/ 1 0 0  ,  T  i  mm ,  290 ,  T empm 
5 0  A X E S  . 001 f 1 0 0 , 0 , 3 © 0 

55  UNCLIP 
60  C SIZE  3 
70  LDIR  0 
o  u  L  U  R  G  b 

SO  FOR  I ~ . 002  TQ  T i mm- , 00 1  STEP  . Q02 

100  MOVE  1,270 

110  LABEL  USING  " HD " ; I  * \ A00 

1 20  NEXT  I 

121  LDIR  90 
130  LORG  5 

149  F 0 R  1=4 O 0  T 0  T e mp m  S T E P  200 
145  MOVE  -.O30*Timm,I 

150  LABEL  USING  " U p B D 11 ;  J 
170  NEXT  I 

1 8  0  L  0  R  G  5 

1 9  0  M  0  V  E  -  ,  O  6  *  T  i  m  rn ,  \  3  O  8  +  T  e  m  p  m )  /  2 

2 0 0  C  SIZE  4 

201  LDIR  90 

210  LABEL  USING  " 20A TEMPERATURE  (BEG.  K>" 

220  LDIR  0 
2  3  O  H  0  V  E  T  i  m  m  /  2  ,  2  3  0 

248  LABEL  USING  "  1  9R"  ;  "TIME  (MILLISECONDS)  11 

2  5  O  M  0  V  E  T  i  m  m  /  2 ,  T  e  m  p  m  +  5  O 

268  C SIZE  6 ,  . 5 

270  LDIR  0 

290  LORG  5 

295  PEN  1 

30O  LABEL  USING  ,,30flM ;  Label  $ 

310  PEN  3 
320  FRAME 
330  PEN  4 
340  MOVE  0,360 
355  OVERLAP 

3 61  D  I  M  Ac  (9 9  )  ,  A  1  (  9 3  ) 

U7U  D  I M  I  st  he  (  7  >  ,  Th (  39,4,7),  Dr<  99,7),  I  st  a<  7)  ,  I  $■%  o  (  7  >  ,  J p  <  7  >  ,  p  t  <  7  )  ,  B s 
, Ak  ( 4 0 ) , Ho  1 d  <  7  > , Kno  <  7  > 

371  DIM  Ah (4) 

372  MAT  READ  Ah 

373  DATA  . 25, 1 , . 5, . 333333333333 

374  Ipropt=l 

3 8  0  READ  T  h  t ,  E  o  p  r ,  L  h  e  1  p ,  T  m  a  x ,  \  \ o  r  d  f 
381  MAT  Hoi d=<Tht ) 

390  DATA  3O0, 1,4,0, i 
40O  ASSIGN  #8  TO  "DATA" 

405  BUFFER  #8 

410  ASSIGN  #7  TO  " DATA 1 " 

411  Eof =6 

420  READ  #7 ; Rhoi , wi , Ti , Du, Tempi , Ada, Xm, D, X i nt 

421  Xi=Xint 

430  INPUT  "GUN  TEMPERATURE? " , Tht 
448  MAT  Th  =  <  Tht ) 

510  MAT  READ  I st a, I st o , Jp , I st he 


C  7> , A 1 ph ( 4  0 ) 


37 


520  DRTR  1,5, 7,8, 9, 9, 9,7, 11, *13, 15, 15 
530  RERD  Ere  ad ,  Er ,  ErO ,  Enew ,  Eos 
540  DRTR  1 , 1 , 1 , u, 1 
550  E  r h  =  u 


"560  H2-1E50 
570  Erh=0 


, 15, 15,4 


4, 4, 4,4, 4,4,  1,5, 7, 8, 9 


9 


580  R E  ft  D  S , 1 r  d  s , N  r d s , K  o  s h 

599  DRTR  @,1,1,0 

600  fl  1  c  p  =  -  4  s  <14.8  +  4*LG  T  ( D  >  >  A  2 
610  INPUT  11 RLCP? 11 ,  R 1  cp 

620  INPUT  "HINT,  TINT?”,  Hint,  Tint 

630  INPUT  “ NRDS, NLfl , IPROPT, LHELP?" , Nrds, HI  a, Ipropt,Lhelp 

640  T=Q=Tcal=0  f 

730  L  i  ni  =  Ln*i  i  n=  1 

740  Tint*B*iO0  * 

750  N  l  a-7 

760  MR T  RERD  Rk,R1ph 

770  DRTR  0,0, 0,0,0, . 10$, . 10S, . 107, . 106, . 105, . 103, . 1@1 , .098, .096, .095, .683, • fl 
.073, .073, .068, .062, .053, .@6, .063, .089, . 107, . 122, . 130, . 1^5, . 140, . 145, . 145 
780  DRTR  0,0, 0,0, 6,0, 0,0 

790  DRTR  0, 0,0, O, u,  . 12t , .  122, . 1  15, . 109, . 103, .096, .09, .084, .078, .072, .066, .061, 

055, .049,  .046, .045,  .@48, . @56f .07, .039, . 113, .  123,  .  135, . 148,  .  145,  . 150, .  150 

800  DRTR  0,0, 0,0, 0,0, 0,0 

320  Rlphfi»=0 

838  FOR  k  1  =  1  TO  40 

840  fil phwsMRX(fll phfrt, R1 ph<K 1 > ) 

850  NEXT  K 1 

860  Tdesir*T1nt 

870  Rlce-422. 5*R1 cp 

880  Hspace=SQR<Du*ftl  phm*3) 

890  Const  stDu/<Hspac  e^Hspac  e  > 

908  FOP  K 1 = 1  TO  40 

910  R 1 ph<K 1 > ^Const  *fi 1 ph C K 1 > 

920  Rk  < K 1 )= Const *Rk <kl > 

930  NEXT  K1 

94G  R1  c  p*fi  1  c  e*Hspa»:  e*Const 

950  PRINT  "THE  SPRCE  I NTERVRL  IS",Hspac€ 

960  INPUT  "INPUT  INSIDE  DIRM,  OUTSIDE  B I RM ” , R i , Ro 
970  Ri=2.54*Ri 

980  Ro=2 . 54*Ro 
990  K  o  u  =  8 - N 1  a 


1 0 O 0  S  h  =  H  u  /  H  s p  ac  e 
1016  D  t  <.  K  o  u  >  =  D  u  *  1 0  0  @ 


1020  Ds (Kou ) =Hspace 
1038  FOR  Kol-Kou  TO  6 
1040  Dt CKol  +  l >=4*Dt  <Kol > 

1 050  Ds < Ko 1  +  1 >  =2*Ds  C Ko 1 ) 

1060  NEXT  Ko 1 

1070  I s t  o < 7  >  =  I  NT  <  < Ro-R i >/Ds<7>  +  .5> 
1030  REM  — CIRCULRR  CORRECTIONS 
1090  FOR  I = K o u  TO  7 

1100  I S  t  Op  = I S  t  O  <1  ) 

1UG  FOR  I  i=l  TO  1st  op 
1120  R  =  Ri +Ds< I >*I i 

1130  DrCIi ,I>=Ds<I)/<2*R> 

1131  NEXT  Ii 

1132  NEXT  I 
1146  Eoe=  1 


1145  R 1 ce=fl Ice/ 12 
1150  FOR  K 7  = 1  TO  4 
1160  L=7 

1170  Kno< 7>=J  =  K7 
1180  GOSUB  Heat  er- 
1190  IF  Eoe  THEN  J7 
1200  FOR  K 6  = 1  TO  4 
1210  L=6 


1220  K  n  o  (.  6  >  =  J  =  K  6 
123  U  G  O  S  U  B  H  e  a  t  e  r 
1240  IF  Eoe  THEN  J6 
1250  FOR  K5 - 1  TO  4 
1260  L  =  5 

1 2 7  0  K  n  o  <  5  >  -  J  =  K  5 
1280  GOSUB  Heater- 
1290  IF  Eoe  THEN  J5 
1300  FOR  K4=l  TO  4 
1310  L~4 

1320  Kno< 4 ) -J -K4 
1330  GOSUB  Heater- 
1340  IF  Eoe  THEN  J4 
1350  FOR  K3--1  TO  4 
1360  L=3 

1370  Kno <3 > =J=K3 
1380  GO S U B  H e at e r 
1398  IF  Eoe  THEN  J3 
1 4 0 0  F 0 R  K 2=1  TO  4 
1410  L=2 

1420  Kno < 2 >  = J  =  K2 
143  0  G  0 o U  £  H  e a t e  r 

1449  IF  Eoe  THEN  J2 

1450  FOR  K 1 = 1  TO  4 
1460  L= 1 

1478  Kno < 1 > = J=K 1 
1 4  8  0  G  0  S  U  B  H  e  a  t  e  r 
1490  NEXT  K 1 
150O  J2S  NEXT  K2 
1510  J3 : NEXT  K 3 
1528  J4 : NEXT  K4 
1530  J5 : NEXT  K5 
1540  J6: NEXT  K6 
155*0  J  7 :  NEXT  K7 
1560  GOTO  1150 

1570  Heater: IF  Kosh< >0  THEN  2690 

1598  DEF  FNDe  1  <  I a>  =Dr  <  la,  < Rk  < I s 1 > -Rk  < I s3 >  > / < 4*Rk  <  I *2  >  >  ; 

1  600  Ii EF  FNF < fi ,  B ,0=01  ( J< > * <R* <  1  -Be )-B~B+C*<  1  tBe >  )  +B 
1610  DEF  FNG  <  fi , B , C , E  > =E  *  F  NF  <  R , B , D  >  +  < 1 -E  > 

1615  E  o p r =  < I p r  o p t  =  1  )  0 R  < I p r o  p t  =  2  >  B  N  D  < J = l > 

1620  Jp<L>=J 
1  6  3  0  R  M  =  0 

1640  Is t art  =  I  s t a  <  L) 

1650  Is t o p = I s t o < L > 

1660  Il  =  Istar-t  +  l 

167G  12= I st op- 1 

1680  Jhe= 1  +  <  J  +  2>  MOD  4 

1690  IF  J= 1  THEN  Ho  1 d < L > =Th < I st op , 4 , O 

1700  FOR  1  =  I  start  TO  I  stop 

171G  L  1  =  T  h  <  I  , J h e , L> / 5  Q 

172G  F*  1  =  L  1  MOD  1 

1725  LI =1  NT (LI ) 

1  730  R1(I)=fll  ph  <  L 1  >  -F  1  *  <  fi  1  ph  ( L  O  -  R 1  ph  <  L 1  + 1  >  > 

1740  He  <  I  >  =  B k  < L  1  >-Pl  *<fik  <L1  >-Bk<L1+l  >  > 

1750  NEXT  I 

1760  IF  L  =  7  THEN  Outside- 
1770  REM  — OUTER  STEP  OF  LRYE*  R 
1780  I s= I st  a(L+ 1 > 

179Q  LI =  Th< Is, Jp<L+l > , L+l >/50 
1  8  0  0  K  =  I  s  t  o  p 

1810  Be-Dr <K, L>-<flc (I st op-2 > -Rk  1 >-Pl >-fik <  L 1 + 1 > >  > / <  4  *  fi  c <Istop)> 

1 820  Th  <  I  s  t  op ,  J  ,  L  >  =  F  N  G  <  Th  <  I  st  op-2 , 4 ,  L  >  }  Th<  I  st  op  ,4,0,  Th<  Is,  Jp<  L  + 1  )  ,  L+  1  >  ,  fih  <  J  > 

1830  Out  side: FOR  K=I 1  TO  12 

1340  De  =  Dr <K, L) - ( Rc <K-1 >-fic <K+i  > <K> > 

1  856  Th  <  K  ,  J  ,  L  >=FNF  <  Th  <  K- 1 ,  .Jhe ,  L  >  ,  Th  <  K ,  Jhe ,  L  )  ,  Th  <  K+l  ,  Jhe  ,L>> 

I860  NEXT  K  39 


1870  IF  L=7  THEN  Th < I st op , J , 7) = <2*Th< I st op- 1 , J, ? > -Th< I st op-2, J , 7 > /2  + 1 . 8 ) / 1 . 566 
1886  REM  --THIS  ASSUMES  K= . 1 1 , H= . 00067, AND  AMBIENT  TEMP  IS  300  K 
1890  IF  L=Lmi n  THEN  A40 
19O0  I s= I st  o<L- 1 ) 

1910  U=ThOs-l,  l,L-l>/50 
1920  P 1 =L 1  MOD  1 

1930  De=Dr<  I  start ,  L>-<Ak  tLl >-Pl  *<ftk  0.1  >-flk  <L  1  + 1  >  )-Ac  <  I  start +  1  >  >/<4*Ac  <Istart>> 
1940  K= I start 

1 950  Th  <  I  s  t  art. ,  J ,  L  )  =FNF  ( Ho  1  d<  L- 1)  ,  Th<  1st  art ,  Jhe  ,L),  Th  <  I  st  art  + 1 ,  Jhe  ,  L  >  > 

I960  IF  Eopr  THEN  PRINT  USING  1955; S, Rv, 6, L>  Th< I st  art , J, L) , Th< 1st  art +  1 ,  J,  L ) ,  Th<  I 
s  t-  ar  t+2,J,L>,Th<Istart+3,J,L) 

1965  IMAGE  BDD2. D, 8DZ, 6D2. DD, 4D, 7DZ. B, 7DZ. D, 7DZ. D, 7D2. D 

1970  GOTO  Helper  _  .... 

1980  A40: S=S+Dt <L) 

1990  IF  Eread  THEN  A58  .  " 

2000  IF  Eos  AND  <Th<  1 ,  J,  LX  1 . 5*Th<5,  J,  L)  )  THEN  A120 
2010  A45:  ThC  1 ,  J,  L>-<2*Th<2,  J,  L>-Th<3,  J,  L>/2>/l  .  5 
2020  Enew=<S>=Tdesi r) 

2.030  IF  Eopr  THEN  PRINT  USING  1965; S, Rv, Sq, L,  Th<  Istart , J , L> , Th< 1st  art  + 1 , J, L) , Th< 
Istart+2, J, L>, Th< Istart+3, J,L> 

2040  GOTO  Helper 

2045  A50:  M-H+l 

2046  READ  #8 ; V , X , Rho , Tgas 

2047  X=X+Xi 

2050  IF  V  > 1 E  4  9  THEN  A70 
2060  A60: Vgas=0 

2070  IF  X>Xi nt  THEN  Vgas»V«X i nt /X 
cd O 8 0  HI  =A  1  c e*Rho*Vgas 
2090  Z=H1*H2/<H1+H2> 

2100  Coef=Z*Sh/flc < 1 > 

2110  Th<  1 ,  J,  L)«<2*Th<2,  J,  L)-Th<3,  J,  L>/2+Coef  *Tgas>/<  1 . 5+Coef  > 

2111  Ry  =  Z* < Tgas-Th < 1 , J, L  )  > 

2112  IF  Eng  THEN  2120 

2113  IF  S>T i mm*1000  THEN  Eng-1 

2119  DRAW  S/1000, Th <1  start , j, L> 

2120  Q =Q+<Ry+R h o v > * D t <L)/2000 

2130  IF  Eopr  THEN  PRINT  USING  1965; S, Rv , Q, L, Th< I st art , J , L > , Th< 1st  art  + 1 , J ,  L )  ,  Th<  I 


£•  t  ar  t  +  2 ,  J ,  L  >  ,  T  h  (1st  ar  t  + 

3, J,L> 

2140 

R  h  o  v  —  R  v 

2150 

GOTO  H  e  1  p  e  r 

2160 

R70:IF  Er  THEN  R86 

2170 

Er ead=0 

2180 

GOTO  ft  6  O 

2190 

ft80 : E  r = G 

2200 

REM - STEP  SIZE 

INCREASE 

2210 

Sh=2*Sh 

2220 

ftlO0:S=S-DtCL) 

2230 

Lrn  i  n=L  i  ni=Mrri=L+ 1 

2240 

I i =Ista<Mm)-l 

2250 

•J.j  =  K  n  o  <  M  m  ) 

2260 

FOR  1=1  TO  Ii 

2270 

Th  < I , J  j ,  Mrn  >  =  T  h  (  2*  I 

-1,4, Mm- 1 j 

2275 

NEXT  I 

2280 

I  s  t  a  (  M  rri )  =  1 

2290 

Kosh=3 

2300 

Kount =0 

2310 

GOTO  2360 

2320 

ft  120:  IF  J<  >1  THEN 

H  45 

2330 

Kount = Kount + 1 

2340 

EO'S=  <L<Lhe  1  p> 

2350 

ON  Kount  GOTO  R45, 

R45 , ft  1 00 

2360 

He  1  per : Eoe=  <  Lm i n  >= 

L  > 

2370 

I  st  art  =  1st  a<  L  > 

2380 

Istop=Isto<L> 

2390 

T  nt  a  x  =  M  ft  X  (  T  rri  ax,  T  h  (  1 

,  J ,  1  >  > 
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2408  IF  Enew=0  THEN  RETURN 

24  1  8  P’R  I  NT  "  Round  no"  ;  I  rds ;  “  T  I  HE=  M 1 888  ;  "  Seconds " ;  LIH<  i  >  ;  "  N  a:*;  t  m p "  ;  I  NT  <  T max  > 
"  Bore  surf  ace  temper at  ure"  ;  Th<  1,1,  U  Pi  e  1  p>  ;  "Tot  al  heat  i  nput 11 ;  .  1  * INT<  1 0*Q) 

2428  Q  =  0 
2430  Enew=0 

2440  IF  Nrds.<  =  Irds  THEN  CftLL  Wai  t  <Nr*ds  •  Tdes  i  r ,  T  i  nt ,  Ret  urn  > 

2450  IF  Return  THEN  RETURN 

24 EG  Tmax=Q 

2478  E o s  =  0 

243  0  I  r  d  s  -- 1  r  d  s  +  1 

2490  E r e a  d = E r = E r g = T r u e 

2580  T  d  e  s  i  r=Td*si  r+T  i  nt 

2510  Sh=Sh/*2 

2520  Lh=,L-l 

2530  FOR  1=1  TO  Lh 

2  5  4  0  J  p  <  I  >  =  4 

2558  NEXT  I 

2560  I=J 

2578  Mmrn  =  L 

2575  E30 : Lh=M- 1 

2  5  y  G  I  s  t  a  C  M  fi’i  rn )  =  1  s  t  h  e  M  fii  rn  ) 

2  598  I  s  =  I  s  t  a  <  M  rn  rn  > 

2608  FOR  Ka= 1  TO  Is  * 

2618  ThC2*K-l ,  4,  Lh)=Th<2*K,  1 ,  Uh>  =  <Th<K,  I ,  M  >+Th<K+l ,  l ,  M> ) / 2 
262  8  I s t o p =  I  s t o < L  h ) 

2630  Ho  1 d <  Lh ) =Th ( I s t op , 4 , Lh  > 

2648  Mmrn=Lh 
2650  1=4 

266G  IF  Lh>=2  THEN  E30 
2678  Lrn  i  n=L  i  rn=  1 
2688  RETURN 
2  6  9 8  K o s h = K o s  h- 1 
2700  E o e  = 1 
271G  RETURN 

2  7  2  O  S  U  B  W  a  i  t N  r  d  s ,  T  d  e  s  i  r ,  T  i  n  t ,  R  e  t  u  r  n  > 

2  7  3  0  R  E  ft  Ii  N  a  d  d  ?  T  a  d  d 
2740  BBT  ft  0,0 

2750  IF  Nadd=G  THEN  BIG  ; 

2760  N rds=N add+Nrds 
2770  Tint=Tadd 
2780  R  e  t u  r  n  =  0 
2790  SUBEXIT 

2800  B 1 0 : IF  Tadd=0  THEN  STOP 
28 10  T des i r=Tdes i r+T add 
2820  Ret urn= 1 
2838  SUB END 
2848  END 
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